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We extend recent work on the resonance structure for post-inflation reheating, providing an 
analytic treatment for models in which both the inflaton and the fields into which it decays are 
massless. Solutions are derived which are valid for either a spatially-flat or spatially-open metric. 
Closed-form solutions are given for the characteristic exponent, which measures the rate of particle 
production during preheating. It is demonstrated that in certain ranges of parameter space, the 
maximum values of the characteristic exponent in an open universe are several times greater than 
the maximum values in a spatially-flat universe. It is further demonstrated that the solutions found 
here by means of a simple algebraic construction match the two previously-known exact solutions, 
which were derived in terms of special functions. 
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I. INTRODUCTION 



In this paper we extend recent work by Greene, Kofman, Linde, and Starobinsky on the resonance structure for 
models of post-inflation reheating in which both the decaying inflaton and the fields coupled to it are massless. 
Their work treats the case of massless fields in a spatially-flat background spacetime. Yet there is increasing evidence 
to suggest that we live in an open universe, and the recent models of open inflation present an interesting way to 
reconcile an open universe with an early phase of cosmological inflation. |2j The formalism of cannot be used 
\f} ■ directly to study the preheating dynamics following a phase of open inflation. Here we extend the work of (!]] to treat 
f*"*- \ preheating in an open universe, and find both qualitative and quantitative changes in the spectra of produced quanta 
between spatially-open and spatially-flat models. In particular, when the initial value of the oscillating inflaton field 
is of the same order of magnitude as the associated curvature scale, there will be no resonant amplification of any 
modes - a result with no correlate in the corresponding spectra for preheating in a spatially-flat universe. Further- 
more, characteristic exponents can reach larger values in an open universe than in a spatially-flat one, signalling a 
Oh, quicker rate of resonant particle production. Because of the conformal invariance of the massless theory, the extension 
provided here also allows study of preheating in Minkowski spacetime with massive fields, with or without explicit 
[ symmetry-breaking terms in the potential. With these extensions, the elegant and powerful methods introduced in 
• may be applied to a much wider class of interesting models. 

The case of massless fields in an expanding spacetime admits a self-consistent analytic treatment for all three of 
• i-h , the time-dependent quantities involved: the oscillating inflaton, the resonantly-amplified decay-product field, and the 
scale-factor for the background spacetime. Whereas previous authors have found closed-form solutions for such pre- 
heating systems in terms of special functions for two particular ratios of the relevant couplings , the authors of Q| 
produced a method for finding solutions for a far broader class of couplings by means of a much more simple, algebraic 
construction. In addition to extending their work to treat preheating in an open universe, we also demonstrate here 
that the Floquet index, which measures the rate of resonant production of decay-product quanta, can often be written 
in closed form, rather than in terms of the "auxilliary functions" found in [|lj , which are defined as definite integrals 
and evaluated numerically. Finally, we demonstrate explicitly that the solutions found by the algebraic construction 
of pj], with only minor modifications, match the exact results found previously in terms of special functions. 

As is by now well known, linear second-order differential equations with purely periodic coefficients have solutions 
which obey Floquet's theorem: the solutions behave as X(t + 2ui) — X(t) exp(iFt), where 2w is the period of the 
coefficients in the governing equation of motion, and F is the Floquet index (also known as the characteristic expo- 
nent). While the form of F will depend on model parameters, it will be independent of t. These solutions can always 
be rewritten in the form X(t) = P(t) exp(/i K i), where P(t) is a periodic function with period 2u>, and /i K is simply 
related to F (see, e.g., j^]). Obviously, when /u, K has non-zero real parts, the solutions develop exponential instabilities. 
Such solutions are said to lie within "resonance bands" or "instability regions," the dependence of which on model 
parameters can often be written explicitly based on an analysis of the form of Floquet's theorem lies at the heart 
of the study of preheating with massless fields, since, for early times after the inflaton has begun oscillating around 
the minimum of its potential, the equations of motion for any fields coupled to it assume the form of second-order 
differential equations with periodic coefficients. It is therefore imperative to develop efficient and accurate methods 
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for evaluating fj, K for a broad class of interesting models. This is the main goal of this paper. 

In Section II, we present the models to be studied and their dynamics. We also discuss the question of appropriate 
vacuum states with respect to which the preheating production of quanta should be measured. Section III presents 
the extension of the method of JlJ for the new cases, and includes explicit solutions of the equations of motion valid 
for early times (before back- reaction becomes significant). In Section IV we demonstrate how to re- write the fi K - 
parameters in closed form for many cases, which facilitates both their numerical evaluation and their comparison with 
known solutions. Appendix A includes an explicit demonstration of the equivalence of the solutions X(t) found in 
Section III with solutions found by other means. In Appendix B it is shown that the /Lt K -parameters calculated in 
Section IV likewise match the exact solutions found previously. 



II. DYNAMICS OF THE MODEL 

The effective Lagrangian density we will study involves two fields: an inflaton <j> and a decay-product field x, both 
of which are massless and have minimal couplings to the Ricci curvature scalar. For the field content of the model, 
we may then write: 



4 , y i2,,2 



(1) 



The spacetime near the beginning of the preheating epoch will assume the form of an FRW metric; for the usual 
models of inflation, this will be a spatially- flat metric (with K — 0), while for the newer models of open inflation, this 
will be a spatially-open metric (with K = —1): 

ds 2 = -dt 2 + a 2 (t) li, ,il.r'd.r ! . 
h^dx^x 3 = dx 2 (K = 0), 

= dr 2 + sinh 2 r (d0 2 + sin 2 9d<j) 2 ) {K = -1). (2) 

For either case, the equations of motion for the fields become 

<f> + ZH4> + Xcf> 3 + g 2 X 2 <t> = 
X + 3Hx-^-L 2 X + .9Vx = 0, (3) 

where dots denote d/dt, H = a/a, and L 2 is the comoving spatial Laplacian operator. We have assumed that the 
slow-roll field (f> is spatially homogenous; that is, we will ignore the spectrum of inflaton fluctuations for now. As 
demonstrated in jlj , the preheating production of inflaton quanta from a massless inflaton with A</> 4 coupling is anoma- 
lously inefficient, as compared with the preheating production of a distinct (massless) boson field. The spectrum of 
inflaton fluctuations may also be studied in this formalism by associating \ with <5</>, and assigning g 2 = A in the large 
N limit of the O(N) approximation, or g 2 = 3A in the Hartree approximation. (See, e.g., 

We may now promote the field x to a Heisenberg operator. Details of this expansion, and of the canonical quanti- 
zation of the resulting operators for both the cases K = and K = — 1, may be found in Q]. For the K = — 1 case, 
we will only track subcurvature modes here; the possibility of amplifying supercurvature modes during preheating in 
an open universe will be treated elsewhere. || The quantum \ field may then be expressed as a sum over modes and 
associated creation and annihilation operators, with 

X(t, r,VL)=J dp, [x P £m(t, r, SI) d plm + H.c] . (4) 

Here djl is the appropriate measure for the integral, including relevant sums over I and m (the specific form depends 
on whether AT = 0orAT = — I, sec Q), a p £ m is a canonical time- independent annihilation operator, and "H.c." denotes 
a Hermitian conjugate. 

In order to study the non-linear dynamics, we will make use of a Hartree factorization. Because we will be studying 
the evolution of the fields before any back-reaction becomes significant, the details of the resonance structures for 
the model of equation ([l ) will in fact remain independent of the specific approximation scheme invoked to study the 
non-linearities (though details concerning how long the preheating phase will last do depend on the choice of such 
schemes). In our case, the Hartree approximation entails replacing the g 2 X 2 term in the equation of motion for <f> 
with its vacuum expectation value, <? 2 (x 2 )- This quantity gives a measure of the growth of the back-reaction upon 
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the oscillating cf> field due to the resonant shift of energy into x modes. Because of spatial translation invariance, this 
quantity can only depend upon time. We will hence study the coupled system for early times, when this back-reaction 
term may be neglected relative to the tree-level terms. 

The equations of motion become much more simple if we rewrite them in terms of conformal time, drj — a~ 1 dt, and 
in terms of the re-scaled fields: 



^ ^ y/Xa(r]) 
£(77) = pf 2 ^))- 



> Xptm(t, r, n) 



i 



01(77) 



Xp(r))Y p i m (r, fl), 



(•5) 



The quantity g 2 Y*(r]) measures the growth of the back-reaction upon <p(rj), and the spatial harmonics Y p i m obey: 



L 2 r p , m = - {f - k) y ? 



pirn ■ 



(6) 



In this paper, we will only be concerned with the continuous spectrum of modes in both the K = and K = —1 cases, 
with the eigenvalue p 2 in the range < p 2 < 00. With these definitions, the equations of motion may be rewritten: 



d 2 _ 1 d 2 a 2 2 S 
drj 2 a drj 2 



1 d 2 a 



9 2 .2 



drf a dr/ 2 ^ 1 



0, 
0. 



(7) 



It is convenient to define the frequencies W p (r)) and uj p (r]) as 



W 2 ( V )=p 2 



K 



W 2 (v) + 9 % 2 {r]) 



1 d 2 a(r]) 
a (77) drj 2 

f 
X 



In this way, W p is the (time-dependent) "natural frequency" for the X p modes in the absence of interactions, and u p 
is the (time-dependent) "natural frequency" for these modes when subject to the g 2 x 2< t )2 coupling. 

One further simplification may be made for the case of preheating with massless fields. As demonstrated in [0, the 
scale factor 0(77), averaged over a period of the inflaton field's oscillations, behaves as if the universe were radiation- 
dominated (even though all the fields are far from thermal equilibrium). Thus, to a good approximation, we have the 
relation 



1 d 2 a{r,) 
a(rj) dr] 2 



K 



(9) 



for the entire preheating epoch. (We thus ignore for now the possibility of additional resonance effects arising from 



the oscillation of a(t); see 



.) This means that the Ricci curvature scalar vanishes during preheating: 



R( V ) 



G 



a 2 (77) [a dr] 



K 



0. 



(10) 



Models of preheating in an expanding universe with massless fields are thus conformal to models of preheating in 
Minkowski spacetime, even though the fields are minimally-coupled to R rather than conformally-coupled to R. 

It remains to consider appropriate initial conditions for the modes ^(77); this is crucial to any study of preheating, 
since the initial conditions for X p come from canonically quantizing the field \ 011 a proper Cauchy surface and 
choosing an appropriate vacuum state around which the quantized field may be expanded. As discussed in [Q, when 
treating only the continuum of modes with p 2 > 0, the proper initial conditions for ^,(77) take the same form for 
both K = and K = —1. There are two independent concerns for choosing an appropriate vacuum state: the 
nonequilibrium interactions amongst the fields (relevant even for discussions of preheating in Minkowski spacetime), 
and the usual ambiguity regarding physical vacua in general-relativistic settings. We choose to study preheating as 
measured against the "adiabatic" vacuum: both the transition from de Sitter spacetime and the "turning on" of the 
interaction between the fields are assumed to occur adiabatically. In the absence of interactions, the transition from 
de Sitter spacetime to the more general FRW spacetime at the beginning of preheating would yield the following 
initial condition for the X p modes at the onset of preheating (taken to be the time 77 = 770): A p (?7o) = [2M / p (t7o)] _1 / 2 
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and (dXp/drj) v — Va — — i[W p (rjo)/2} 1 ^ 2 . These modes would represent the free-particle "adiabatic" states ||. If the 
interactions were also turned on adiabatically beginning some time before 770, then these initial conditions would be 
replaced by: 



1 



V /2 ^p(?7o) ' V d V 



dX r . 



(11) 



These give the initial conditions for the "adiabatic" particle states for the nonequilibrium dynamics at the onset of 
preheating, at r\ — rj . With this choice of vacuum state, the number of produced quanta per mode may similarly be 
written in the same form for both K = and K = — 1 01: 



N„ 



ujpiri) 



P 



dX v 



dq 



(12) 



Only when the \ field is properly quantized will this yield the number of "adiabatic" -state quanta produced per mode 
relative to the initial Fock space vacuum. It is clear that when p lies in a resonance band and X p ~ e^" 11 , then 
N p ~ e 2 ' tK,J , giving rise to the very efficient transfer of energy from the oscillating ip field into x quanta. 

Finally, there is the issue of renormalization. With our choice of vacuum (and hence of initial conditions for X p ), 
the term (X 2 (r]o)} is formally quadratically divergent. This ultraviolet divergence may be treated in either of two 
ways: we may undertake a formal renormalization, as done, for example, in || (see also |fTo| ), or we may use the 
"physical" criterion that the energy density contained in the vacuum fluctuations at the beginning of preheating 
must be less than that contained in the classical inflaton field, <f>, as done in (J]. In practice, this amounts to setting 
g 2 T,(r]o) <C (Pq(t]o). Calculationally, the divergence will not be a problem, because preheating generically produces 
many quanta in low-p resonance bands, and we will only be 'following' the evolution of these exponentially-growing 
unstable modes, so that in practice we will not have to evaluate a full integral over < p < 00. We may now turn to 
the evolution of the fields during preheating. 



III. EVOLUTION OF THE FIELDS DURING PREHEATING 

Using equation (^|), the equation of motion for jp in equation (Q) may be solved for early times, when g 2 E(r/) may 
be neglected. Setting 770 = 0, the solution is ||[|,[l]]: 

(p(ji) = cn LfpJ+K n , </V v 2 ^o + i o) , (13) 

where cn(u, v) is the Jacobian cosine function, and this solution is valid for tp^ > 2\K\, appropriate for chaotic inflation 
initial conditions. Note that for non-zero K, the opposite limit, ipo <C \K\, will not produce any resonant preheating 
effects: in this case, (p{rj) reduces to the non-periodic form ^2(77) ~ ipoCOsh(\K\i]). 

The solution in equation (|l3|) also describes the early-time evolution of a massive inflaton for models in Minkowski 
spacetime: in this case, K < corresponds to preheating with an explicit symmetry-breaking potential, while K > 
corresponds to preheating with a massive inflaton but without symmetry breaking. If we scale the fields by the 
inflaton mass, then these reduce to K = ±1. || For the remainder of this paper, however, we will restrict attention 
to preheating in an expanding universe with massless fields, and hence K will equal either or —1 depending on the 
spatial curvature of the metric. 

At this point, we could follow the procedure of jUJ] and substitute this solution for </?(?/) into the equation of 
motion for the modes X p (r]). After much calculation, explicit forms for the modes may then be found in terms of 
special functions for certain values of the ratio of couplings g 2 /X- Instead, we will adopt the methods developed in 
jlj to produce exact (early-time) solutions for X p by means of algebraic construction. This approach may be applied 
whenever the couplings satisfy <? 2 /A = n(n + l)/2, with n a positive integer. With this ratio of the couplings, the 
equation of motion for X p in equation (Q) takes the form of a Lame equation of order n. 

The authors of study the case of K — 0. With the help of the following definitions, we may expand their method 
to the case of non-zero K: 

l=\l<Pl+ K > U = 1V , z = cn 2 (u,i/), (14) 
where v = ipo/^2(ipQ + K) = ipo/ v^T 2 - The Jacobian elliptic functions obey the relations (see, e.g., pi]] ) 
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sn 2 (u, v) + cn 2 (u, v) = 1 , dn 2 (u, v) + v 2 sn 2 (u, v) — 1 , — cn(u, v) = — sn(u, v) dn(u, v) 

du 



(15) 



so that 



du 



-2[(l-v 2 )z + {2v 2 -1) 



Dz 2 



Given the bound <~p\ > 2\K\, the modulus v obeys 



- < i/ < 1. 
2 ~ 



(16) 



(17) 



Deviations of v 2 from 1/2 indicate effects from the non-zero spatial curvature; the limit v 2 — ► 1/2 is the same as 
the limit ip§ 3> \K\, that is, the limit when the initial amplitude of the oscillating inflaton is much greater than the 
associated curvature scale. 

The Jacobian cosine function cn(u, v) is periodic with period AK{v) (as measured in 'u' units of time), where K(y) 
is the complete elliptic integral of the first kind. This means that z(u) (oc (p 2 (u)) has period 2K(v) = 2u>, and hence 
that Floquet's theorem applies to the equation of motion for X p in equation (0). In terms of z, the equation of motion 



for X p becomes: 



Af{z)X'; + 2f{z)X' p + (V + 2^ 2 ^ X p 
where primes denote d/dz, and we have defined 

2 _ 2/2 
K =p /7 . 

Finally, we have defined the function 



/(z) = (1 - u 2 )z + {2v 2 - l)z 2 



v 2 z\ 



(18) 



(19) 



(20) 



With these definitions, we may study equation ( fl8| ) for X p for any value of the couplings g 2 /X = n(n+ l)/2, following 
the general approach of JlJ . 

Because the differential equation for X p is second-order, there will exist two linearly-independent solutions, which 
we can label Ui and Ui (suppressing the index p for the moment). Consider the product M(z) = Ui(z)U~2(z). After 
straightforward algebra, and making use of equation (p"8|), one can show that this function obeys the third-order 
differential equation: 



2f{z)M"'{z) + 3f'(z)M"(z) + 



f"{z) + 2[K 2 + 2v 2y -z 



M\z) + 2v 29 —M(z) =0. 
A 



(21) 



When g 2 /X — n(n + l)/2, this equation may be solved in terms of polynomials in z of order n, which we will label 

M (n) {z): 



(n) Jn-i) 



(22) 



and we will set = 1 for all n. Because z(u) is periodic in u with period 2uj = 2K(v), M/ n \(z) will also be periodic 
in u. We will give some explicit examples below for n = 1 and n — 2. With these solutions for Mr n \(z) in hand, one 
can then find solutions for Ui(z) and ^2(2), using the equation for X p and the Wronskian relation between U\ and 
U 2 - 

In general, for a differential equation of the form 



d 

dx 



A(x) 



dy_ 

dx 



B{x)y = 0, 



the Wronskian for the two linearly-independent solutions will be proportional to l/A(x). In our case, we have 



d_ 

dz 



dz 



B{z)X p = 0, 



(23) 



(24) 
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and hence 

U[U 2 - UiU' % = 4===, (25) 

v/O) 

where the constant C(„) will be a function of k and v and will depend on the order n. Using this Wronskian and the 
fact that Ui(z)U~2(z) — M(z), we may write solutions for the two mode functions: 

Cf n i f dz 

yW)M (n) {z) 



U lA z) = ^^^exp ± ^> / I . (26) 



We have inserted a normalization factor N, which is set to unity in [Q. Instead, we will set 

N =\M (n) (z = l)\- 1 / 2 , (27) 

since z = 1 at u = 0. With this normalization, when the solutions are written as functions of u, they satisfy 
£^i,2 = 0) = 1; and, as demonstrated in Appendix A for the case n = 1, they match the exact solutions found 
previously, by very different means. 

With Ui^{u = 0) = 1, the normalized mode solutions U\ 2 do not obey the proper initial conditions at u = r\ = 0, 
equation ([□]) . We must instead take linear combinations of these solutions p,|| : 



X p (u) ' 



2^2^(0) 



jU[{u = 0)7 V jU{{u = 0) 



(28) 



where primes here denote d/du. Unless such linear combinations of Ut t % are taken, the number operator N p in equation 
( |l2| ) will not measure particle-production with respect to the physically-relevant initial vacuum state. Numerically, 
however, if one makes the approximation of only following the exponentially-growing unstable modes, then N p will 
only shift by a constant as one shifts the normalization of Ux^ and the linear combination for X p . 

The solutions will be exponentially unstable whenever the exponent in equation ( p6| ) has a non-zero real 
part; we will use this to determine the bounds on the resonance bands in terms of n and v. Furthermore, given 
the periodicity of <f(u) and the quasi-periodicity of Ui^{u), the authors of jjj demonstrate that the characteristic 
exponent may be written as a definite integral: 



= _C(n) f 1 dz 

2a, ./o yW)M (n) {z) 



We will present means for further evaluating fi K given this expression in Section IV; in particular, we will provide 
a transparent demonstration that the definite integral in equation ( p9[ ) is always purely real, so that all resonance 
structure comes from the behavior of C^ n y We will demonstrate explicitly in Appendix B that this form for fi K 
does indeed match the exact solutions found earlier in |3|Q] for the two specific cases already studied. We will 
show in Appendix A that the solutions in equation ( |26| ) obey Floquet's theorem, and can be written as f/i^w) = 
P(±u) exp(=F/x K w), with P(u + 2uj) — P(u). 

By plugging the solutions for in equation (^6|) back into the equation of motion for X p , one finds (dropping 
the subscript (n) for the moment): 

C 2 = M' 2 f - 2M"Mf - M'Mf - k 2 M 2 - v 2 n(n + l)M 2 z. (30) 

It is straightforward to confirm that all of the z-dependence in this expression cancels exactly for any order n, leaving 
the simpler relation: 

Cl^-K^a^y-a^Hl-v*). (31) 

Note that equations (|l|), @), @, (|2|L_and (||) all reduce to the forms in (T) when K = (y 2 = 1/2). 
For g 2 = A (that is, n = 1), equation fl21|) may be solved, yielding: 

k 2 1 
M w (z)=z--z + — -2, 

v A C 2 (1) = -k 6 + 2k 4 (1 - 2^ 2 ) - k 2 (1 - bu 2 + hv A ) - v 2 (\ - 3is 2 + 2v A ). (32) 
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When K = (v 2 = 1/2), these become M (1) (z) -> z - 2k 2 and C 2 1} -> k 2 (1 - 4k 4 ). Thus, in this limit, there will 
exist one single resonance band for positive k 2 , determined by when C%\ > 0, and given by 

< k 2 < -. (33) 

This matches the resonance structure found in both P,|lJ. It is interesting that in the opposite limit, with non-zero 
K and ip\ ~ (y 2 — > 1), C 2 ^ — > — k 6 — 2k 4 — k , and there will not exist any resonance bands for k 2 > 0. In 

fact, in an open universe, there will exist a minimum value of ipo for any given ratio g 2 /A below which no modes with 
k 2 > will be resonantly amplified. This issue will be treated further in || . 

For the next-simplest case, g 2 — 3 A (n = 2), it will be easier to write the solutions in terms of the coefficients of 
the expansion in equation ( p2[ ) (recalling that = 1): 

M (2) (z) = z 2 + af z + 4 2) , 

4 2) = -^[« 2 -4(l-2^)], 

a ( 2 ) = JL [ K 4 _ 5/t 2 (1 _ 2 ^ 2) + 4 _ 25iy 2 + 25l/ 4] ) 

^ 2 ). (34) 
k 4 ) , (35) 



(36) 

matching the resonance structure found in Note that in the opposite limit, v 2 — > 1, 

C ( 2 )--^ 2 (« 4 + 4) 2 , (37) 

again revealing no resonance bands for k 2 > 0. 

We now turn to the treatment of the characteristic exponents, 



c ( 2 2) = - K 2 (4 2) ) 2 - a ^d- 

When K — 0, this reduces to 




revealing the presence of a single resonance band for positive k 2 : 

\ < k 2 < V3, 



IV. EVALUATING THE CHARACTERISTIC EXPONENTS 

As noted above, the authors of |Q found a form for the characteristic exponent, fi K : 

c^r^dz 

2 W 7o V7(I)A/ ( „ ) (z) 

In this section, we will rewrite the definite integral in a way which demonstrates explicitly that the integral is real for 
all n and v 1 so that the resonance bands are determined entirely by where the constant C7 n ) has non-zero real parts. 
In addition, the integral may be evaluated in closed form in many cases. Such closed-form solutions for fi K may then 
be evaluated without numerical integration, and will facilitate comparison with known solutions, as discussed in the 
appendices. 

With the substitution 

z = l- sin 2 9, (39) 



the indefinite integral in equation (Gq) becomes: 



=2f r , * — (40) 



V?(ijM (n) (z) J [l-^ sin 2 0]V2 M(n)(sin 2 
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Thus, given 1/2 < v 2 < 1, both the integral over dz and the integral over d9 will be manifestly real for all n and v. 

The form of the integral over d9 in equation ( [40|) further encourages comparison with the incomplete elliptic integral 
of the third kind (see, e.g., p2|): 

d9 [1 - n sin 2 0] ~ [l - sin 2 a sin 2 0] ~ 1/2 . (41) 

Whenever the function M(„)(sin 2 9) may be factored into n real roots, the integral over dz may be written explicitly 
as a sum over n distinct II- functions. In the cases when M(„)(sin 2 9) cannot be so factored, then the (principal value 
of the) integral over d9 in equation (|40| ) will still be well-defined over the limits of integration (z : 0, 1) — > (9 : n/2, 0). 
Consider the case when M(„) (sin 2 !?) can be factored into n real roots; writing x = sin 2 9, we have: 



where the n constant coefficients may be determined by the set of n equations: 

n I n 

Ea lift 

A ( / 9 2 - 1 +/3 3 " 1 + - + ^ 1 )+c.p. = 0, 
r>! (/3 2 - 1 /3 3 - 1 + /^Vr 1 + - + C-i^n 1 ) + C -P- = o, 

Di (fc 1 ^ 1 x ... x + c.p. = 0, (43) 

where c.p. denotes all cyclic permutations. Note that when the constant portion of M( n )(sin 2 f?) (that is, a^) is 
non-zero, all of the roots Pi will be non-zero, and hence their inverses (3^ will always be well defined. In general, 
both the constants Di and Pi will be functions of n, k 2 , and v. We will give explicit examples of the factorization of 
equation (^2|) for n = 1 and n = 2 below. 

In the cases when M(„) (sin 2 9) admits such a factorization, then we may write the indefinite integral in the exponent 
of U\.2 as: 

f dz " 

~ / m — ? n j / \ = 2 E DiU (Pi 1 ' arcsin ^ 1 ~ z \ arcsin v ) ■ ( 44 ) 

J y/f[z)M {n) {z) ^ 

Using the fact that II(rt; \a) = and H(n; tt/2 \a) — H(n\a) (the complete elliptic integral of the third kind), we 
may further evaluate the characteristic exponent, /i K : 

Mk = -^M f 1 dz = -^yAn(/3- 1 \arcsinz,). (45) 

For cases in which a given (3J 1 > 1, one may always rewrite II(/3~ 1 \ arcsin v) in terms of a sum of complete elliptic 
integrals, each of which is purely real (see |l2|). This closed- form solution for /i K , valid whenever M( n )(sin 2 9) has the 
n non-zero, real roots Pi, may thus be evaluated in terms of well-known functions, without any need for numerical 
integration. 

Consider the simplest case, n — 1. In this case, M^^sin 2 9) may always be factored as needed: 

M (1) (sin 2 9) = - sin 2 - ^ + J_(i _ v 2 ), (46) 



or, 



Pl x =D X = - - -. (47) 

1 — V — K 



This yields 



8 



^ = fl T jVfl T 2 -\arcsin^. (48) 

Equation ( p8[ ) reveals that preheating in an open universe can be more efficient than in a spatially-flat one. The 
maximum for (i K when K = (^ 2 = 1/2) is /i K = 0.147 at k 2 = 0.228 (matching the result found in Q). Yet /i K grows 
in the range 1/2 < v 2 < 0.76, reaching a maximum of fi K = 0.224 at k 2 = and j/ 2 = 0.76. This qualitative difference 
in the spectra will be treated further in || . We will demonstrate below in Appendix B that the analytic solution for 
fi K in equation (^) matches the exact solution found previously for n — 1, with v 2 — 1/2. 
For n = 2, we have (again using x = sin 2 6): 



M {2) 



(x) =x 2 -x(2 + o[ 2) ) + (l + af 5 + 4 2) ) , (49) 



or (dropping the superscript (2)) 

A.2 - ^ 
^1,2 = I 



2 + ai T \/(ai) 2 -4a 5 



2 (1 + oi + a 2 ) \/(ai) 2 - 4a 2 



±2± ai + v/( ai ) 2 -4a 2 . (50) 



Note that when (ai) 2 < 4a 2 , the definite integral in fi K will still be purely real; it simply will not be expressible as a 
sum of complete elliptic integrals. In the limit v 2 — * 1/2, however, the condition (ai) 2 > 4a 2 is equivalent to k 4 < 3 
(see equation ([54j)). And in this limit, the single resonance band extends only to k 4 < 3, so for the entire resonance 
band, the roots and coefficients in equation ( |50| ) will all be purely real. 

As with the n = 1 case, preheating in an open universe with n — 2 can be more efficient than in a flat one. In 
this case, the enhancement of the resonance can be far more dramatic: whereas fi K reaches a maximum value when 
K = (y 2 = 1/2) of [i K = 0.036 (at k 2 = 1.615), fi K rises sharply over the range 1/2 < v 2 < 0.90. At v 2 = 0.72, fi K 
reaches a maximum of 0.097 (at n 2 = 0.64), while at v 2 — 0.90, fi K reaches a maximum of 0.193 (at k 2 = 0). This 
yields a characteristic exponent over five times greater than the maximum reached in a spatially-flat universe. We will 
demonstrate analytically in Appendix B that for n = 2 and in the limit v 2 — > 1/2, the closed- form solution obtained 



from equations (M5j) and (50) again matches the exact solution. 



V. CONCLUSIONS 

The methods developed in |l| for the analytic study of the resonance structure of preheating with massless fields are 
powerful and highly efficient. We have extended their work here to cover both models with non-zero spatial curvature 
(relevant to studies of open inflation) , and models in Minkowski spacetime with or without explicit symmetry-breaking. 
By considering a physically well-motivated choice of initial vacuum state, we have also fixed the overall normalization 
of the mode functions, which provides a simple form for the number of x quanta produced per mode during preheating. 
With this choice of normalization, the solutions for the mode functions found here by simple algebraic construction 
match the solutions found previously in terms of special functions (see Appendix A). 

Beyond the relative ease with which solutions for the mode functions may be found using the methods of [EJ, 
there is an added benefit of this approach over that taken in, e.g., HQ: the characteristic exponent, /i K , which 
determines the rate at which particles are resonantly produced during preheating (with \nN p ~ 2fj, K u), may be 
evaluated independently of finding the full solutions for the mode functions, f/i, 2 (-z). We have shown here that in 
many cases of interest, these /i K parameters may in fact be solved for exactly in closed form, in terms of well-known 
functions. Written in this form, the characteristic exponents, calculated directly from equation (pOf), match the 
solutions found previously (see Appendix B). 

Of course, this method, though simple and efficient, is still limited to the cases in which the ratio of the couplings 
satisfies g 2 /X = n(n + l)/2, with n a positive integer. The authors of (j]] show, by means of numerical integration 
for arbitrary positive g 2 /A, that /i K is not a monotonic function of the couplings; for slightly different values of g 2 /X, 
the resonance can be much stronger than it is with some of the integer-valued ratios studied here. Still, it is useful 
in general to have analytic solutions in hand, especially if one wants to compare different models, such as K = 
versus K — — 1, or, for Minkowski spacetime, > versus m 2 < 0. Such analytic comparisons have revealed here, 
for example, that when K = — 1 and the initial amplitude of the inflaton's oscillations falls below a minimum value, 
no modes with k 2 > will be resonantly amplified, even though there is no associated threshold in the K = case. 
Moreover, when K = — 1, certain regions of parameter space yield rates of resonant particle production over five times 
greater than the corresponding rates in a K = universe. These topics will be treated further in 
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APPENDIX A 



In this appendix, we demonstrate that the solutions Ui^[z(u)] in equation (p6|), with the normalization of equation 
(p7|), match the exact solutions found in terms of special functions. We will examine the simplest case, n — 1, which 
was studied in 0J. In terms of the parametrization of equation (|T4|), and setting v 2 = 1/2, we may follow the same 
steps as in [|| to arrive at the solution: 



(51) 



and Ui,u(u) = U K (±u). Here $i(x) are the Jacobian theta functions, ui = K(v) is the complete elliptic integral of the 
first kind, and Z(x) is the Jacobian zeta function (see, e.g., |ll]Ji2]| ). The quantity v is related to k 2 by 0: 

p(2ujv + uj') = -k 2 , (52) 

where p(x) is the doubly-periodic Weierstrass function. We have also used u>' = iK'(v), where K'(v) = K(v'), the 
complete elliptic integral for the complementary modulus, v' = (1 — v 2 ) 1 / 2 . |pd| , ^2| We may rewrite the v-dependence 
in terms of Jacobian elliptic functions, using 



sn 2 (u, v) 



v 2 sn 1 (u + lo' ,v) 



= 2p(u + uj') + 1, 



where the last expression holds for v 2 — 1/2. Making these substitutions, equation ( |52|) may be rewritten: 

1 - 2k 2 = sii 2 (2luv, v). 



(53) 



(54) 



Inside the resonance band, with < k 2 < 1/2, we thus have < v < 1/2. 

In order to demonstrate the equivalence between the solutions in equations (|2^) and (|5l]) , we will first examine the 
term outside of the exponent in equation (vw: 



J\M {1) (z)\ = Vz - 2k 2 = [sn 2 (2ww, v) - sn 2 (w, v)] 1/2 . 
Again using equation (^3|), the following relations between the theta functions [p"i"| 



i?i ( x + — t ) = i exp 



-ITT [ X + —T 
4 



1?4 | .r + — r J = / (\\p 



(where r = lu'/lu), and the relation (for v 2 = 1/2) 



—vk x H — r 
4 



p{u) 



1 

AJ 2 



gj(0) giM) 



-1 2 



we may rewrite equation (^5|) in terms of theta functions: 

1 i?i(0) 1 



\M {1) {z)\ 



1/2 



The normalization of equation (27) may be similarly re-expressed: 

N = 1 



Vl - 2k 2 sn(2wu, i/) ' 



(55) 



(56) 



(57) 



(58) 



(59) 
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Finally, using the relation between the squares of the theta functions J13| : 

- M^Kiv) = Mv + z)Mv ~ ^l(o), 

we find 

j? 4 (0) r. / u \ . t u mV- 



Now we may turn to the exp onent in equation (E™ . 



1?! 



2oj) V 2lo) 



(60) 



(61) 



Using equations (|4J), (|47|), (|53|), and (|54), with z/ 2 = 1/2, we may write the exponent of equation (|26j) as 



C(i) 



= — C(i)Diil I -sn 2 (2uw + u/); arcsin V 1 — i \ arcsin(^) 



Because < isn (2u>v + u>') < i/ , we may rewrite the incomplete elliptic integral as follows |lz 



II(n; y>\a) = 5 a 



1 { Mv + P) \ , 



09) 



05) 



(62) 



(63) 



with 



e = arcsin(n/sin 2 a) 1/2 , /3 = F(e\a)/2w, 

y = F(ip\a)/2u} , S 1 = [n(l - u) _1 (sin 2 a ~ nT 1 ] 1/2 



(64) 



where F((p\a) is the incomplete elliptic integral of the first kind. Here we have used the V conventions for the 
arguments of the theta functions of [0] . After some straightforward algebra, equation (|62]) may then be rewritten: 



(i) 



dz 



u & 4 (v) 
2uj i? 4 (w)' 



Vf{z)M {1) (z) 

after use has been made of equation ( |56"| ) and the fact that -d^—x) = i? 4 (x). Finally, using |ll| 



2ujZ(2ujx), 



(65) 



(66) 



and combining equations (61) and (|65|), we find that the solution for Ui[z(u)] found in Section III, equation (|26|), may 
be rewritten as 



Mv) Mjh) 



(67) 



Given the periodicity of the theta functions, this is now in the form: Ui(u) = P(u) exp [— fJ, K u], with P(u + 2uj) = P(u). 
The second independent solution can then be written ^(u) = P(—u) exp [/x K u], and represents the growing solution. 

This completes our demonstration that the solutions (for n = 1) found here by means of very simple algebraic 
construction, following the method of match the exact solutions found earlier by much more difficult means ||. 
The comparison for n = 2 between the solutions found in Section III and those found earlier in Q| follows very 
similarly. 



APPENDIX B 



Whereas in Appendix A we demonstrated that the solutions 1/1,2(2) for the mode functions match the earlier known 
solutions, in this appendix we will demonstrate that the crucial quantities, fi K , as calculated directly from equation 
( [i"5| ) also agree with earlier known solutions. The difference between these two demonstrations lies in the fact that the 
solutions 1/1,2(2) are defined in terms of an indefinite integral over z, as in equation (|26j); the characteristic exponents 
fji K , however, are defined as definite integrals. As noted in the conclusions, one benefit (among many) of following the 
methods of |l| is that the fx K parameters may be evaluated independently of the mode functions ^1,2(2), which is not 
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possible if proceeding as in EM . 

Consider the case n = 1 with v 1 — > 1/2. Equation ( fl8|) then becomes: 



Written in this form, (3 1 1 > 1, so we may rewrite the elliptic integral as follows i jjjfl : 

n(/3f V) = u) - IL(N\a) = -8iLoZ{e\a), (69) 

with N = /3i sin 2 a, and (Si and e defined as in equation (|64|), now in terms of iV instead of n. Then, using equation 
we have 



e = arcsin \J\ — 2k 2 = arcsin (sn(2cjw)) . (70) 
Noting that when tp = arcsin (sn(x, v)) p2| , 

Z(tp\ arcsin(j/)) = Z{x\v 2 ) = Z(x), (71) 

we find that 

fi K = Z(2ujv) (72) 

for n — 1 and v 2 — > 1/2. This matches the result found in 

For n = 2 and v 2 — > 1/2, we may compare with the solution found in 

u -JLf^ + flfs)^ (73) 
where the constants a and 6 are defined by the relations 



1 2 1 / 1 „ 

<p(a = k \ 1 « , 

8 v; 6 2 V 3 ' 



p{b) = -\^ + \^l-\^. (74) 
Writing these in terms of Jacobian elliptic functions, using the identities employed in Appendix A, we have 



sn 2 (a - uj') = 1 - -k 2 - Wl - -k 4 , 
V ' 3 V 3 ' 



sn 2 (fo-c') = l-^ 2 + ^l-^4. (75) 



Examination of equations (|34|) and ( |50| ) in the limit i/ — > 1/2 reveals that 

z^r 1 = —2?-" — a = ^ sn2 («)' 

sn z (a — a/ J 2 



Then 



3 fK 2 + l 



An 2 V k 2 - § 



-1 + -k 2 - W 1 - -rt 4 II -sn 2 (a)\ arcsin^) + 1 - -k 2 - \ 1 - ~k 4 II -sn 2 (&)\ arcsin(^) 



(77) 
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Since these j3 i 1 each satisfy < (3 i 1 < 1/2, we may again use equation (|(33|) for each of the two complete elliptic 
integrals. In each case, the \n.{§i(y + /3)/i?4(?/ — /3))-term vanishes, leaving, after much algebra, 

(78) 

Thus, the closed- form solution for found in Section IV, equation (|45|), matches the solution found previously in Q| 
for n = 2. 
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